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We examine the exchange Hamiltonian for magnetic adatoms in graphene with localized inner 
shell states. On symmetry grounds, we predict the existence of a class of orbitals that lead to 
a distinct class of quantum critical points in graphene, where the Kondo temperature scales as 
Tk oc \ J — Jc\^^^ near the critical coupling Jc , and the local spin is effectively screened by a super- 
ohmic bath. For this class, the RKKY interaction decays spatially with a fast power law ~ 1/FC . 
Away from half filling, we show that the exchange coupling in graphene can be controlled across 
the quantum critical region by gating. We propose that the vicinity of the Kondo quantum critical 
point can be directly accessed with scanning tunneling probes and gating. 

PACS numbers: 71.27.+a,73.20.Hb,75.30.Hx 
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Graphene is a single atomic sheet of carbon atoms with 
elementary electronic quasiparticles that behave as mass- 
less Dirac fermionsfl'|. The Kondo effect has been re- 
cently observed in graphene 0, and the formation of 
a Kondo screening cloud around a magnetic adatom is 
quantum critical at half filling0,[^, crossing over at weak 
coupling to the standard Fermi liquid case, when the 
DOS is locally restored by disorder or gating effects[3]- 
The Kondo resonance in graphene is also strongly sen- 
sitive to the position of the adatom in the honeycomb 
lattice, where the interplay of orbital and spin degrees of 
freedom may give rise to an SU(4) Kondo effect [8|. 

In this letter, after establishing a generic one-level ex- 
change interaction Hamiltonian for magnetic adatoms in 
graphene, we show there is a symmetry class of orbitals 
in which quantum interference between the different hy- 
bridization paths leads to a fixed point where the Kondo 
temperature Tk cx: \ J — JcY ^ scales with the mean field 
exponent v = 1/3, with J as the Kondo coupling near 
criticality. In the v ^ l/i class, graphene behaves as a 
super- ohmic bath for the local spin and the RKKY in- 
teraction is strongly suppressed, decaying spatially with 
a fast power law ^ \/EJ . Furthermore, we show that 
the exchange coupling in graphene can be controlled by 
gating. This effect opens the possibility of exploring the 
proximity to the Kondo quantum critical point (QCP) in 
graphene directly with scanning tunneling probe (STM) 
measurements i9r;12j . 

We start from the graphene Hamiltonian, T-Lg = 



].{'Ri)ba-{'R.j) -f h.c, where a,b are fermionic 



operators on sublattices A and B, respectively, t ~ 2.8 
eV is the nearest neighbors hopping energy and a ^fj, 
labels the spin. In the momentum space. 



ng = -t^ 0p4,p^(T,p + h.c, 



(1) 



where (j)p = J2^=i e*^ '*', and ai — x, ai2 — —x/2+V5y/2, 
and as — —x/2 — ^/3y/2 are the lattice nearest neighbor 



vectors. 

In the presence of a localized level, the problem is de- 
scribed by the single impurity Anderson Hamiltonian 13 



IJ, H=ng+Hf+Hu+Hv, where = J2a eo^^/,< 



IS 



the Hamiltonian of the localized electrons, with fif^^ — 
f^fa as the number operator, and eo is the energy of 
the local state measured relative to the Dirac point, 
Jiif — Ufif^-^fif_^ gives the electronic repulsion in the lo- 
calized level, and Hy describes the hybridization between 
the local level and the graphene electrons. 

The adatoms in graphene can sit for instance on top 
of a carbon atom, where the hybridization Hamilto- 
nian is T-Ly = V^^a\,{Q)fc^ -\- h.c, or in the hollow 
site in the center of the honeycomb hexagon, where 
the adatom hybridizes with the two sublattices, T-Ly = 

E.Eti [K,.4(a.) + H,z6t(-a,)] fa + h.c, with F,,, 
{x = a, b) representing the hybridization strength of the 
localized orbital with each of the three surrounding car- 
bon atoms sitting on a given sublattice. In momentum 
representation 



(2) 



per 



where Vap = V , and Vb.p = for a top carbon site, 
say on sublattice A (A-site). When the adatom sits in 




Figure 1: Representation of (a) an s-wave and (b) an /-wave 
orbital, when the adatom sits in the center of the graphene 
honeycomb hexagon, (c) Substitutional impurity in a single 
vacancy. In the three cases, the adatom hybridizes equally 
with the neighboring carbons on the same sublattice. 
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the center of the hexagon (ff-site), or else for a sub- 



stitutional impurity in a single vacancy 15| (S'-site) , the 
hybridization function is Vx.p = K^.^e^^ '*'. On H- 

sites, for an s-wave orbital, Vxa = V, giving V^.p = V(f>p, 
whereas for an in-plane /-wave orbital, as shown in 
Fig. lb, where the orbital is odd in the two sublattices, 
Vfo,,; = -Va.i = V, resulting in -Va,p = Vb^p = F0p. In 
the case of an s or in-plane /-wave orbital on an S'-site 
on sublattice A, Va,p = 0, and Vf,_p = V(j)p, whereas for 
a substitutional impurity on a _B-site, Va^p = V(l)p, and 
Vb,p = (see Fig. Ic). 

Diagonalizing the non-interacting part of the Hamilto- 
nian Ji in the A, B sublattices, 



Tig - 



E 



(7,p*-'Q(,cr,p : 



(3) 



where Ea{p) = ai|0p| is the graphene tight-binding spec- 
trum, with a = ± labeling the conduction and valence 
bands, and c±^cr.k = (fe(T,k±</'k/l'/'k|acr,k)/'\/2 are the new 
operators in the diagonal basis. The hybridization term 
in the rotated basis is 



-± p,(T 



where 



(H,p+aK:p0;/i0pi)/(V2v^). 



(4) 



(5) 



In particular, O^p — l/\/2 when the adatom is on top 
of an A-site, 6f p = a(j)p/ {^/2\(j)p\) on a B-site, and 

6a,p [<f>P + (-l)''a'/'pVl'/'pl]/V2 when the adatom 
sits on an H-site, where 7 = for an s-wave orbital 
and 7 = 1 in-plane /-wave orbital. In the substitutional 
case, Qafp = 't'p/V^ for an impurity on sublattice A, and 
ef^p = a(/i*2/(|0p|V2) on sublattice B. 

For all possible symmetries, the orbitals of adatoms 
sitting on 5* or sites can be classified among those 
that either break or preserve the Ca^ point group sym- 
metry of the triangular sublattice in graphene. Since \(j)p\ 
scales with \uj\/t, the orbital level broadening, A(a;) = 
Ea.p l©a,pl'^(^ - ai|<^pl) is either A(w) ex V^p{lj) 
for orbitals that explicitly break the C^^ point group sym- 
metry, in which case |6q,pI scales to a constant near 
the Dirac points, where p cx |a;| is the graphene den- 
sity of states (DOS), or else A{uj) cx V'^ p{uj)\uj\'^ /t'^ , 
for C'sv invariant orbitals, when |Oc(,p| oc \(j)p\ scales to 
zero at small energy. The first class of orbitals, where 
A(w) cx |a;| (say, type I), represents the standard case 
of ohmic dissipation [31, and is described for instance 
by adatoms on top carbon sites, by Ei{dxz,dyz) and 
E2{dxy,dx'2^y2) representations of d-wave orbitals and 
fxz-^, fyz^, fxyz, fzix^-y^) orbitals in H/S sites. The sec- 
ond class, where A(a;) cx (type II), represents a 
new class of super-ohmic dissipation [l6j, and is described 
by s, dzz, fz3, fx(x^-3y^), and fyCix^-y2) orbitals in H or 



S sites (see Fig.l), where the adatom hybridizes equally 
with the three nearest carbon atoms on a given sublat- 
tice. On physical grounds, this new class emerges from 
quantum mechanical interference between the different 
hybridization paths in the honeycomb lattice, as the elec- 
trons hop in and out of the localized level. As we will 
show, these two classes of orbitals are described by two 
distinct types of Kondo QCP. 

The Anderson Hamiltonian in graphene can be sepa- 
rated in two terms, % = Hq + Hy, and then mapped 
into a spin exchange Hamiltonian through a standard 
canonical transformation, H = e^He'^ — H + [S,T-L] + 
[S,n]] + ... , where 5 = - Ep,a. ^[(1 - ^/,-.)(eo - 

ai|(^p|)~"^ +n/_-<j(eo + U - at\(f>p\)~^]<da,pci,pcrfa - h.c. , 
which results in a Hamiltonian that is quadratic in V to 
leading order, H = ^0 + [S, Hv]/^ + 0{V^)17]. At large 
the exchange Hamiltonian is given by 

■^e = - J X! E ®"k0a'k' S • cl,^^,^^,aCa^^^i, , (6) 

where a = (ai, (72, fXs) are Pauh matrices and 

Jip)^V^U/[{eo-^i){eo + U-^l)]<0, (7) 

is the exchange coupling defined at the Fermi level, p,. 

The validity of the exchange Hamiltonian ^ is con- 
trolled by the ratio A(eo)/|eo — /i| ^ 1, when the valence 
of the localized level is unitary (and hence, the local spin 
is a good quantum number) and perturbation theory is 
well defined in the original Anderson parameters [isj . In 
graphene, where A(a;) cx ttV"^ p{uj){\uj\/t)^ , (ry = 0, or 2), 
with p{oj) = \uj\/D'^, and D ^ 7eV as the bandwidth, this 
criterion becomes | J| - V^/{p - Eq) < DH"^ /{TT\eo\^+'i). 
When the level is exactly at the Dirac point, eo = 0, the 
level broadening is zerofl^ and the exchange coupling 
|J| ~ V^/p has no upper bound and can be shifted by 
gating towards the strong coupling limit of the Kondo 
problem, J — >■ oo, when the Fermi level is tuned to the 
Dirac point, p — >■ O"*"!!^. Since the experimentally ac- 
cessible range of gate voltage for graphene on a 300 nm 
thick Si02 substrate is p & [—0.3, 0.3] eV, the exchange 
coupling of a magnetic Co adatom, for instance, with 
V = 1 eV and eg = —0.4 eV, can be tuned continu- 
ously in the range between | J| g 1.4 — 10 eV. This effect, 
which is allowed by the low DOS in graphene, brings the 
unprecedented experimental possibility of controlling the 
exchange coupling and switching magnetic adatoms be- 
tween different Kondo coupling regimes in the proximity 
of a QCP, as we show in Fig. 2a. 

Since the determinant of the exchange coupling ma- 
trix in Eq. ([S]), det[jQ,Q'], is identically zero, the 
exchange Hamiltonian ^ can be diagonalized into a 
new basis where one of the channels decouples from 
the bathji^l. The eigenvalues in the new basis are 
-'^u.k.k' = Jj^a^a-k^ak' and Jy = 0, and hence, the 
generic one-level exchange Hamiltonian ^ maps into 



the problem of a single channel Kondo Hamiltonian, 
T-Le = — 2 X^k k' "^«,kk'S • Sk,k', where s is the itinerant 
spin, regardless the implicit valley degeneracy. 

In the one-level problem, the renormalization of the 
constant J due to the coupling of the local spin with the 
bath is given by: J' = J - 2NsJ^piD){D/t)'^ 5D/D, af- 
ter integrating out the high energy modes with energy D 
at the bottom of the band, where Ns — 1, 2 describes the 
number of sublattices the adatom effectively hybridizes. 
Since a DOS in the form p{uj) oc has a scaling dimen- 
sion r, where r = 1 in graphene, the restoration of the 
cut-off in the "poor man's scaling" analysis requires an 
additional rescaling J' [{D + SD)/DY+''iJ'[A], which 
results in the beta function 



/3(J) 



dJ 
dlnL> 



-ir + 7^)J-2NsJ^piD)iD/t)^. (8) 
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Figure 2: (color online) a) Kondo coupling vs. gate p for 
U/t = 1 and V/t = 0.3. The dot illustrates a typical value 
for the critical coupling Jc- b) Kondo temperature Tx/t vs. 
J for Csv invariant orbitals. Red (light) curve: p. = 0; black: 
fi/D =0.05, 0.1, 0.15, 0.2, 0.25, and 0.3. The line J ~ Jc sets 
the crossover scale between the Kondo weak coupling (WC) 
and strong coupling (SC) regimes. Inset: Tk /t vs. J/Jc — 1 
near the QCP {p = 0), in log scale. Red (light) line: type II 
orbitals {v = 1/3); black: type 1 — 1) (see text). 



The renormalization group (RG) flow leads to a in- 
termediate coupling (IC) fixed point at Jc — —(r + 
■q)t^ /\^Nsp{D)D^], which separates the weak and strong 
coupling sectors. For type I orbitals (ohmic bath), one 
recovers the usual IC fixed point Jc — —r /[2Nsp{D)\^, 
whereas for type II (super-ohmic bath, 77 = 2) Jc ~ 
—3t^/{2NsD) in the Dirac case (r — 1). In graphene, this 
new fixed point describes a one-channel Kondo problem 
in the presence of an effective fermionic bath with DOS 
p oc 1^1'', where f = r + rj ~ 3. Since the tree level 
scaling dimension of the hybridization V in the Ander- 
son model is (1 — r)/2, the case f = 1 corresponds to an 
upper critical scaling dimension, above whichfr > 1) 
is an irrelevant perturbation in the RG sense [22|. In this 
situation, fluctuations are not important near the QCP, 
and the critical exponents are expected to be mean-field 
like, in contrast with the marginal case (r — 1), where 
mean field cannot be trusted[23[. 

The RG analysis derived from the exchange Hamil- 
tonian @ can be verified directly from the hybridiza- 
tion Hamiltonian In the large iV limit near the 
critical regime, singly occupied level states are en- 
forced at the mean field_ level through the constraint 



0[21|, with = 2 in the spin 



1/2 case. The minimization of the energy d{'H)ld\ — 
gives A = J^^dujn{uj)lm[G-ff{u!)T,-ff{uj)], where 

n{uj) = [e("-'^)/^ -I- is the Fermi distribution, T is 
the temperature, and Gffij) = — (r[/CT(r)/J(0)]) is the 
/-electrons Green's function, Gff{uj) — [uj ~ eq — \ + 
^ff{u:) + ^0+]-\ Effiu;) = E.,p |e«,pP G^^M is 
the self-energy of the /-electrons in the presence of the 
graphene bath, where A(a;) = — Iml]//(a;) defines the 
level broadening and G°;^(aj) = (w — a|0p| -I- iO+)~^ is 
the retarded Green's function of the c-electrons in the di- 
agonal basis, G° (r) = — (T[ca(T)cJ,(0)]). In the critical 
regime, where A = /i — eg = y^/| J|, the Kondo tempera- 
ture in graphene is extracted to leading order in V from 



the equation 

1 
J 



N 



2 ^ ai|0p 

p,a ' ^ 



■ tanh 



atl 



2Tk 



(9) 



In the Dirac cone approximation, the Kondo tempera- 
ture for orbitals of type II (77 = 2) is 



Tk = {D/2) (1 - Jc/Jf'^ 



(10) 



at half filling, where Jc = — 3t^/ {NgND) is the same crit- 
ical coupling derived from the RG equation ([S]). Away 
from half filling, Jc defines the crossover between the 
Fermi liquid weak coupling regime, at J/ Jc <C 1, where 
Tk = (|/i|/2)exp[Z?V(3|Mp)(l-Jc/J(M) + 3AiV-D')],and 
the strong couphng regime, for | J| > | Jd w (2/A^s)eV, 
where Tk ^ {D/2)[l - Jc/J{p) + ip^/D'^Y/^, as shown 
in Fig. 2b. At the critical coupling (J ~ Jc), 



TK = {D/2)\SpyD 



2:1/3 



(11) 



and the fingerprint of the QCP at ^ = can be observed 
in the scaling of the Kondo temperature with p in the 
vicinity of the QCP, at J ~ Jc. This scaling can be mea- 
sured in STM, where the signature of the Kondo effect is 
manifested in the form of a Kondo resonance in the DOS 
at the Fermi level, for T <Tk- 

In Fig. 2b, we numerically calculate the scaling of the 
Kondo temperature in tight-binding. For type II orbitals, 
the u = 1/3 exponent in the Kondo temperature, Tk oc 
I J — Jcl", found in the linear cone approximation persists 
above room temperature, up to TK/t 1 (red curves). 
In the more standard ohmic case, for spins on top carbon 
sites (black curve of the inset), the scaling is linear {i/ = 
1) at the mean-field level. 

Tracing the conduction electrons in the exchange 
Hamiltonian the RKKY Hamiltonian of a spin lat- 
tice in graphene is Hrkky = -J^ J2ij Xlf Si -Sj , where 



4 



a) 



DC 
^ -1 



-2 




b) 



o 



Figure 3: (color online) a) x(^) distance R, along a zigzag 
direction (in lattice units), for ^-sites (black triangles), H- 
sites (blue squares), and S sites for spins on the same sublat- 
tice (red circles). Solid lines: fi — 0; dashed: ^ = t. Inset: 
\Xij\ plot in a log scale. Orange (light) guide line: 1/B?\ 
black: I/FC . On the right: x(l) for ^ site spins, plotted in 
the graphene BZ at b) = and c) /j. = t. Red (dark) regions 
represent x{^) > and blue (light) regions x(.^) < 0- 



Rj) is the spin susceptibility, with i,j 
indexing the local spins, and x,y = A, B, H, Sa, Sb la- 
bel the position of the magnetic adatoms in lattice. In 
momentum space. 



X^'^(q) 



j^x,y n[Ea> (p + q)] - n[Ea{p)] 



'a,a',p,q ~ ^ap^ap^a'p+q'^a'p+q- 



(12) 



where 7W^^„ _ . 

For spins on the same sublattice, Ai^^ = 
1/4, whereas on opposite sublattices = 
aof' (/)p0*^q/(4|0p| |(/)p+q|), in agreement with Ref.fSH'l, in 
the Dirac cone limit. For an iJ-site 12611 . 



M 



HH 



\^a,p\ \^a',p+q\ : 



(13) 



where 10 



a,p I — 



[l + (-l)'^aRe((/)3)/|0p|3] for or- 



bitals of type II; for S'-sites, M^^^^ = |0p|2|0p+q|2/4 
for spins on the same sublattice, and Al-^-*-^^ = 
aaVp0p+q/(4|?!'p||0p+q|) for opposite ones. 

In Fig. 3a, we show the spacial decay of the RKKY 
interaction on the lattice for type II orbitals on A, H 
and Sa site spins. At half filling, the RKKY interac- 
tion is always ferromagnetic for same sublattice spins, 
substitutional or not, and antiferro for spins on opposite 
sublattices 2J, |25|. The H case on the other hand, is 
ferromagnetic for nearest neighbor spins and antiferro- 
magetic at longer distances (blue squares) . In the H and 
S cases, the interaction is short ranged and decays with 
a fast power law ~ in contrast to the known l/R^ 

decay in the A site casejij, [1^ [H, 111], as shown in the 
inset of Fig. 3. This fast decay is consistent with the 
case of carbon nanotubes, where the RKKY interaction 
decays with 1/R for top carbon sites and with for 
isotropic orbitals on H sites [29|. 

Fig. 3b and 3c display the magnetic peaks in the sus- 
ceptibility in the A-site case for /i = 0, and fj, = t. For 



IJ. < t, x{q) has a strong ferromagnetic forward scattering 
peak around the center of the BZ {q = 0), and six sub- 
dominant antiferromagnetic peaks at corners of the BZ. 
Exactly at /i = i, a strong peak emerges at the M point 
due to the nesting of the Van-Hove singularities (VHS) 
of the graphene band (see Fig. 3c), where the DOS di- 
verges logarithmically. This peak reverses the ordering 
pattern of the RKKY interaction in comparison to the 
/i = regime in all studied cases, as shown in the dashed 
lines of Fig. 3a. When /j, is at the VHS, the interaction 
between spins on same (opposite) sublattices, substitu- 
tional or not, is always antiferromagnetic (ferro). In the 
same way, the RKKY interaction in the H site case be- 
comes antiferromagnetic for nearest neighbor sites and 
ferromagnetic at long distances. 

In conclusion, we have derived the one-level exchange 
Hamiltonian for magnetic adatoms in graphene and 
shown the existence of two symmetry classes of mag- 
netic orbitals that correspond to distinct classes of Kondo 
QCP. We also showed that the exchange coupling can be 
controlled across the quantum critical region with the 
application of a gate voltage. 
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Erratum: Kondo Quantum Criticality of Magnetic Adatoms in Graphene 
[Phys. Rev. Lett. 106, 016801 (2011)] 

Bruno Uchoa, T. G. Rappoport, and A. H. Castro Neto 

In Ref. [H we have predicted the existence of two different classes of quantum critical points for the Kondo problem 
in graphene, which was shown to correspond effectively to the problem of a localized spin 1/2 coupled to a fermionic 
bath with electronic density of states p{uj) oc , with r = 1 or 3. Here we point out that the mean field exponent 
i/ = 1/3 derived within the slave boson approach for the class of orbitals of type II [f — 3) is incorrect. 

Above the upper critical scaling dimension of the Anderson model (r > 1), the intermediate coupling fixed point is 
non-interacting and describes the level crossing between singlet and doublet states with the trivial exponent z/ = 1 . 
Albeit fluctuations do not play a role in the critical behavior for r > 1, the critical theory is not of the Ginzburg- 
Landau type and the validity of the mean-field slave boson equation of state (9) breaks down in the f — i class, 
invalidating Eq. (10), (11), and the inset of Fig. 2b for the case of type II orbitals. 

All the other results of the paper remain valid, including the spin exchange Hamiltonian in Eq. (6) and the 
prediction of a fast power law for the spatial decay of the RKKY interaction for type II orbitals 1/R'^). 

We note that since hyperscaling is not obeyed for f > 1, the scaling prediction for the Kondo temperature with the 
chemical potential, /i, in the quantum critical region, Tk oc \fi\, can be violated 3]. In the situation where the scaling 
prediction fails, the criticality in the r — 1 and r = 3 classes can be in principle distinguished. That will be verified 
with NRG methods elsewhere. 

We acknowledge M. Vojta for many helpful discussions. 
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